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Abstract 

In the framework of the non relativistic quark model, an exhaustive 
study of radiative transitions in mesons is performed. The emphasis is 
put on several points. Some traditional approximations (long wave length 
limit, non relativistic phase space, dipole approximation for El transi- 
tions, gaussian wave functions) are analyzed in detail and their effects 
commented. A complete treatment using three different types of realistic 
quark-antiquark potential is made. The overall agreement with experi- 
mental data is quite good, but some improvements are suggested. 

1 Introduction 

Quantum Chromodynamics (QCD) is today the only reliable theory for de- 
scribing strong interactions. There exist many systems that can be used as a 
laboratory for exploring and testing the properties of this basic theory. Among 
them, the meson and baryon sectors have deserved a lot of investigations, essen- 
tially because they are very easily produced. However they belong to the non 
perturbative application of QCD and thus are not easy to be described from 
first principles. Despite many improvements in recent years, on both theoretical 
and computational sides, the lattice gauge calculations are still not completely 
reliable and cannot explain the whole bulk of known properties, even for the 
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simplest systems such as the mesons, which consist of a valence quark-antiquark 
pair. 

This explains why so many various phenomenological approaches have been 
developped in order to describe the non perturbative part of QCD. Among 
them the non relativistic quark model (NRQM) has met with an impressive 
number of successes [Q. The puzzling question is that it still works even in 
situations where it is expected to fail; there exist certainly some deep reasons 
for such a behaviour although they have not yet been clarified precisely (see 
[||). Basically the NRQM needs to solve a Schrodinger type equation with 
two-body quark-quark (or quark-antiquark) interactions. In recent years, the 
determination of the interaction between constituent quarks has reached a high 
degree of sophistication and the whole spectra of mesons for instance can be 
accounted for in a rather satisfactory way || . 

However one does know that the description of the spectra is a necessary but 
not a sufficient condition for aiming at a good explanation of non perturbative 
QCD. In particular several very different potentials can give rise to spectra of 
the same quality. One needs other observables in order to test more precisely 
the resulting wave functions. A possibility is the study of static properties, such 
as magnetic moments or charged square radii. But more sensitive observables 
concern the transitions between various states or production mechanisms (which 
rely essentially on the same dynamical operators). One can think for instance on 
meson decays under strong forces (a resonance giving two or several mesons) or 
the decays under electroweak forces (a resonance producing a photon or leptons 
in the final channel). The advantage of this last kind of transition is that the 
transition operator is known perfectly well and thus it is easier to disantangle 
the drawbacks coming from less well known strong interactions through the 
meson wave function. 

In fact this statement is not completely true in the NRQM. Being a phe- 
nomenological theory, NRQM deals with effective degrees of freedom, the con- 
stituent quarks, and a pure Dirac form of the quark-photon vertex in the tran- 
sition operator is questionnable. Moreover, even in the traditional approach 
of radiative transitions (decay of a resonance into a resonance of lower energy 
plus a real photon) several types of approximations are of current use; the effect 
of these approximations can hide the necessity of using a more sophisticated 
vertex for the quark-photon coupling. Most of those approximations originate 
from the formulae widely used in atomic or nuclear physics which are simply 
translated in the meson sector. Let us mention the dipole approximation for El 
transitions, the long wave length approximation (LWLA) and a non relativistic 
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phase space factor. 

Although these approximations are fully justified in atomic or nuclear physics, 
it is not obvious that they continue to work when applied to mesons. Indeed in 
this sector the transition energy is typically E 1 — fc 7 = 0.1-0.5 GeV, while the 
size of the source is roughly R — 0.5-1 fin = 2.5-5 GeV -1 , so that the long wave 
length condition k 7 R << 1 is not really justified. Comparing the photon energy 
to the mass of the emitting meson also convinces us that a non relativistic phase 
space is probably not suited. Moreover the fact that the electrons in an atom 
or the nucleons in a nucleus have the same mass is not true in the case of some 
mesons, and new phenomena can appear. 

During the seventies and eighties, a lot of works has been done on radiative 
transitions for mesons (and also for baryons but we are not so interested in this 
sector here). At the very beginning they were studied in the vector dominance 
model [f|, ||, then the quark model was introduced either in the framework of 
MIT bag model @, @, ||, the non relativistic quark model |], [T(J 0, 2 body 
Dirac equation (jJJ (0| @ or some relativized phenemcnological quark models 
| fl5[ [l(|. But even in the most complete and nice works, as jl5| or jyj, there 
is always an approximation or an inconsistency which plagues the results or 
forbids to draw precise conclusions. Most of relativistic models suffer of a bad 
treatment of the center of mass, relativized model do not treat with equal care 
the quark-quark potential and the electromagnetic operator. Moreover in many 
papers, a non relativistic phase space or a long wave length approximation are 
used and we will see that this is not justified in the meson sector. In addition, 
only very few studies concern the totality of the known experimental data but 
focus on very specific transitions (light quark sector or heavy quark sector or 
even more restricted sets). 

The aim of this paper is essentially twofold. First, within a very precise 
framework namely the NRQM, we want to investigate deeply and with a par- 
ticular care the effects of all these approximations as compared to an exact 
treatment. Second, we compare, on an exhaustive list of transitions and using 
an exact treatment, how different meson wave functions (obtained with different 
potentials giving meson spectra of similar quality) influence the results at the 
level of radiative transitions. We have in mind to see whether it is necessary to 
modify the quark-photon vertex; our study is thus a first necessary step before 
undertaking a more difficult and ambitious program. 

In this paper we will consider all the radiative transitions (which are suf- 
ficiently reliable) that are reported in the particle data group booklet because 
we want an exhaustive analysis. The experimental data can be gathered into 
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several groups : 

• the transitions allowed by LWLA; they are essentially Ml transitions 
( 3 Si -> 1 So or 1 S -> 3 Si) and El transitions ( 3 Pj -> 3 5i or 3 5i -> 
3 Pj); there is also the particular El transition corresponding to the decay 
61 (1235) -» 7T 7 ( l P x -» ^0) 

• the transitions forbidden by LWLA; they are scarce but interesting : they 
correspond to 3 Pj — > 1 5o, 3 Si — > 3 Si and x Pi — > 3 Si transitions. 

The paper is organized as follows. In the next section we present how the 
meson wave functions are obtained and also the different quark-antiquark po- 
tentials that we are studying. In the third chapter we recall the formalism 
necessary for the description of radiative transitions putting the emphasis on 
the general treatment and the differences corresponding to the various approxi- 
mations that we want to discuss. In the fourth chapter our final expressions for 
the total widths are summarized. In the fifth chapter the results are presented 
and the effects of each approximation are analyzed in detail. Conclusions are 
relegated to the last chapter. 

2 Description of mesons 

In the NRQM, the meson is considered as a two particle system : a (constituent) 
quark of mass mi and an antiquark of mass mi submitted to a potential V(r), 
so that the corresponding Schrodinger equation writes : 

P 2 

[mi+m 2 + Y + v ( r )] I *<*) = m <* I **) CO 

where [i is the reduced mass, p the relative momentum, and m a the total mass 
of the resonance. This last quantity, as well as the wave function | *f? a ), of 
course depend on the choice for the potential. The ordinary quarks u and d can 
be considered as an isospin doublet; they are noted generically as n. 

2.1 The potentials 

In this paper we will consider three different types of potential, the so-called 
Bhaduri's (BD) potential [0, AL1 and API potentials |l8|| . For the purpose of 
our analysis, it is not necessary to introduce very sophisticated forms including 
spin-orbit, tensor forces, ...whose effects are not so important and which compli- 
cate quite a lot the formalism. We limit ourselves to potentials containing two 
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different structures : a central term and a hyperfine term; this is the minimum 
requirement to get reliable wave functions : 



V(r)=V c (r) + V h (r)a 1 .<r 2 (2) 

In this expression, we forget about the color term which is always the same 
whatever the meson under consideration. 

The central term V c (r) contains a short range coulomb part (remnant of one 
gluon exchange) and a confining term : 

VAr) = -- + ar p -C (3) 
r 

BD and AL1 potentials exhibit a traditional linear confinement (p = 1) 
while API potential uses p = 2/3, a power best suited to Regge trajectories in 
non relativistic dynamics. In each potential the hyperfine term is short range 
(remnant of the Dirac factor of the Fermi-Breit approximation). For BD it is of 
Yukawa type: 

k e XP (-r/ro) 

TO1TO2 TTq 

For AL1 and API it is of gaussian type: 

= 2* eM-rVfo) (5) 

3mi?Ti2 77 L > /Tq 

but, contrary to BD, the size is mass dependent through 

r = A(2/i)- s (6) 

The parameters for the potentials are gathered in table ([j]). The meson spectra 
obtained with these potentials are rather good and have been presented else- 
where jf8j. It is important to stress that they also give quite good results in 
the baryon sector; so we have some confidence that they describe the quark 
dynamics in a meson in a rather satisfactory way. 



2.2 Meson wave functions 

Because of the rotational invariance, the meson wave function is written as: 

I Vilsj) =»7i(l,2)[* ni (l,2)xs(l,2)]j (7) 

where 77/(1, 2) is the isospin wave function with total isospin /, XsO-i 2) is the 
spin wave function with total spin S and $„l(1, 2) the space wave function with 
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Table 1: Parameters for the three used potentials, defined by All the 

units are in power of GeV. 





BD 


AL1 


API 


m u = m d = m n 


0.337 


0.315 


0.277 


m s 


0.600 


0.577 


0.553 


m c 


1.870 


1.836 


1.819 


nib 


5.259 


5.227 


5.206 


K 


0.520 


0.5069 


0.4242 


a 


0.186 


0.1653 


0.3898 


P 


1 


1 


2/3 


C 


0.9135 


0.8321 


1.1313 


K 




1.8609 


1.8025 


ra 


2.305 






A 




1.6553 


1.5296 


B 




0.2204 


0.3263 



orbital momentum L and radial number n. Spin S and orbital momentum L 
are coupled to total angular momentum J, but are nevertheless good quantum 
numbers. The color wave function is forgotten since it does not play any role in 
the formalism; the flavor content of the meson is also forgotten since the elec- 
tromagnetic operator does not change the quark flavor. The magnetic quantum 
numbers are not indicated here although the projection of the isospin plays a 
role since the operator is a mixing of isoscalar and isovector terms. We will 
come back on this point later on. 

We will use the expression of the space wave function both in position rep- 
resentation and in momentum representation. Thus we write : 

(r|$„ L (l,2)> = R nL (r)Y L {f) 

(p|$„ L (l,2)) = RnL( P )Y L (p) (8) 

We use the same notation for the radial wave function R n L in both represen- 
tations if there is no risk of confusion. In fact the momentum representation is 
very often more convenient to deal with the matrix elements of the transition 
operator. Moreover an approximation of the exact wave function in terms of 
gaussian functions will be particularly well suited to compute easily difficult 
quantities. In this case we write: 
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N 



Rul{p) =P L Y. c * exp(-A,p 2 /2) 



(9) 



i=l 



and a similar expression for the wave function in coordinate representation. In 
general the number of gaussian terms needed in (|^) to achieve convergence is 
rather weak. Just to give an idea of the quality of such an expansion in the case 
of AL1 potential, the masses m w and m p as function of N are presented in Table 
(|)and compared to the exact ones (obtained by solving the differential equation 
with the Numerov algorithm). In the same spirit we plot the corresponding 
wave functions in figure (|l|). Although there is some differences between the 
N=l and N=2 cases, one sees that very rapidly the approximation (||) can be 
identified to the exact solution. In practice we perform our calculations with 
iV=5 and consider the corresponding wave function as the exact one. Let us 
also point out that, since the basis states are not orthogonal, different sets of 
(cj, Ai) parameters can be used with equal success. In fact, we determine these 
coefficients by two different methods : i) minimisation of the state with respect 
to these parameters, ii) best fit of the exact wave function with a form like (^). 
For N—5, both procedures give exactly the same results. Thus in the rest of 
the paper we qualify the wave function given by (|^) with N = 5 as the exact 
wave function. 

3 Radiative transitions 

A number of points are already well known, and we do not want to spend too 
much time on them. We will focus our attention essentially on new aspects or 
on formulations that are discussed in details later on. Everywhere in this paper 
we employ natural units h = c = I. 

3.1 Transition operator 

We start with the non relativistic expression of the electromagnetic transition 
operator between an initial meson state and a final meson state plus a real 
photon of momentum k, energy E=| k | and polarisation e(k, A). We adopt, as 
usual, the Coulomb gauge, and we normalize the plane waves in a quantification 
volume V. 




(10) 
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Table 2: Convergence properties, with three different quark-antiquark poten- 
tials, for the 7r and p masses as function of the number of Gaussian functions 
NG in (Q). This exact value is obtained with the Numerov algorithm. 





BD 


AL1 


API 


7T IG 


0.252050 


0.193702 


0.192776 


2G 


0.157151 


0.140359 


0.140285 


3G 


0.141707 


0.138436 


0.138391 


4G 


0.138835 


0.138220 


0.138192 


5G 


0.138233 


0.138188 


0.138168 


exact 


0.138186 


0.138057 


0.138044 


p IG 


0.780344 


0.771104 


0.772196 


2G 


0.778425 


0.770731 


0.770166 


3G 


0.778402 


0.770033 


0.770156 


4G 


0.778375 


0.770013 


0.770036 


5G 


0.778374 


0.770006 


0.770029 


exact 


0.778347 


0.769958 


0.769988 



2 

M = ^2^- exp(-ik ■ r,){2p l - ia, x k) (11) 

2—1 

The summation runs on the two particles of charge a and mass rrii present in 
the meson. The first term in M is known as the electric term and the second 
one as the magnetic term. 

3.2 Transition amplitude 

The initial meson of mass m a is at rest and has angular momentum J a M a 
(coupling of L a and S a ), isospin I a Mj a . The final meson of mass m& has a 
total momentum K&, angular momentum JbMj, (coupling of Lb and Sb), isospin 
IbMi b . Inserting M of ( |Tl"| ) between the corresponding wave functions (in the 
rest frame of the meson Pi = - P2 = p) in momentum representation (^) (com- 
pleted by the center of mass plane wave), one gets the transition amplitude: 

M A ^ B = 5 Ki< - k [M%i B + M%l B ] (12) 
M%l B = Jd 3 p$ B (p-— ^_k)[2p- J a 1 xk]^^$^(p) (13) 
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M 



(2) 



d 3 p$* B (p 



'" 1 -k)[-2p-i(7 2 xk]i^$,(p) 



mi + ni2 ' ' 2m2 

where the subscripts refer to the particle number, see figure (||). 

We must say a word about the charge and its relation to isospin. One can 
perfectly well ignore isospin degrees of freedom and speak only of flavor wave 
function and its symmetry properties. In this case ej is simply a physical quan- 
tity, the charge corresponding to the flavor of quark i; we understand < > 
appearing in ( |l3| ) as ej. On the other hand, one can also introduce isospin 
degrees of freedom for convenience, in particular, it makes things easier in the 
case of neutral mesons with flavor mixing or mixing angles; for isospin U=0 one 
still has < ei > = ej. However for isospin tj=l/2, the charge is different for 
each member of the multiplet so that it becomes an isospin dependent operator 
ej = ±l/6 + (ti) z (+ for quark, - for antiquark). In this case on must understand 
< ej > as the matrix element of this operator on isospin wave function r\j. If 
we adopt this philosophy, a very simple calculus, based on Racah algebra, gives : 



< e, > 

< ei > 



— e i^I a I b ^MrM Ih 



>M, Mi. 



gfr.* + {-l) 1 '^- 1 ' V3(/« + 1/2) 



<I a M Ia lO | I h M h > 



1 1/2 1/2 

h la lb 



h = 1/2(14) 



< e 2 > = 5 



Mi a Mi, 



~Slah + {-l) h+tl - 1/2 VWa~+W) 

6 



< I a Mj 10 I I b M h > 



1 1/2 1/2 

tl la lb 



t a = 1/2 



3.2.1 Long Wave Length Approximation 

Due to the recoil term in meson B wave function, the space term appearing 
in (^) is not easy to calculate. A widely used approximation is the long wave 
length approximation (LWLA) which consists in putting k=0 in the argument of 
<&b (this is equivalent to replace exp(— ik-r) by 1 in coordinate representation). 
It is just a matter of Racah algebra to disantanglc the spin and space degrees 
of freedom in (|l3|). It is pleasant that, in this case, the electric and magnetic 
part, which are of different parity, cannot couple the same states; this is why we 
speak about electric and magnetic transitions. It is also convenient to calculate 
the covariant spherical components /i = —1, 0, 1 of the vector M. 
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The electric transitions change the parity of the state, but not the spin. One 
has, for the /i component : 



(M^ B ) M = S Sa ,sM b ,-k < J a M a ln | J b M b > F e (A -» B) (15) 
with the expression F e , independent of the magnetic numbers: 

F e (A - B) = (_i)S.+J.+L*+i ( !£L^ _ <^>)Lj a < L a 010 | L b > 

mi m2 




/ d P R nbLb (p)R naLa (p)p 3 (16) 
Jo 



where we have introduced the usual notation J=\/2J + 1. 

The magnetic transitions do not change parity nor orbital momentum. Since 
the operator contains a spin dependent part, the decoupling is a little bit more 
complicated but straightforward anyhow: 



(M^V = 5 La , Lb 5 Kb ^F m (A -> B) 

< Iflff I 1/* >< JaAfoli/ | J 6 Af 6 > (17) 

V,<T 

with the expression -F m , independent of the magnetic numbers: 



mi fni 

{ S a L a J a 1 f 1 1/2 1/2 1 

Ji 1 S 6 J 1 1/2 S Q S b [ 



(18) 



/>OC 

/ dpR nbLb (p)R naLa (p)p 2 
Jo 



Of course, when the gaussian expansion of the wave function (|!j) is employed, 
the dynamical integrals are analytical. 

3.2.2 Beyond long wave length approximation 

The big advantage of using the wave functions expressed on gaussian terms ([}[) 
is that the treatment of the general case can be dealt with rather simply. Since 
taking N=5 is equivalent to treat the exact wave function, the following treat- 
ment solves exactly the problem. In fact an individual term in the expansion is 
of the form $l(p) = exp(— Ap 2 /2)3^i(p) where 3^i(p) = P L Yl{p) is the usual 
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solid harmonic. Just to illustrate the procedure, let us consider only one term in 
the expansion (one A quantity for meson A, one B quantity for meson B) with 
a unit amplitude c A — c B = 1). The argument in meson B wave function is now 
a linear combination of p and k. Such a combination in the solid harmonic can 
be treated using the formulae : 



i 

y lm (a Pl + bp 2 ) = ^C/ 1 a il 6 / - il [^ 1 (p 1 )^_ il (p 2 )] ; 
h=o 



C l = 4tt(2/ + 1)! 

h V (2/i + l)!(2(/-Zi) + l)! 1 ' 



and 



[yiMyiAp)]im = B i hl2 p i ^- i y lm ( P ) 



b ». - { o o ; ) < 2 °> 

The same combination appearing in the exponential is a quadratic form 
which can be diagonalized in order to get rid of the non diagonal terms. Let us 
just summarize our conclusions. 

Here are some auxilliary quantities (i = 1,2 refers to particle number): 

2 ~ (m 1+ m 2 )(A + B) [Zi) 

z {i) = m ^ A . F «) = Dx {i) z d) 

(mi + m 2 )(A + B) 

If the masses of the quark and the antiquark are the same there is no need to 
distinguish the quantities x, z, F and further simplifications arise which we do 
not want to comment. The transition amplitudes in the general case can be 
converted to a more appropriate form: 



M£L B = [d*qeM-Dq 2 -F( 1 h 2 )[y* Lb (c l -zWk) XSb } Jb 

LTYL\ J 

[2q-i<7i x k\[y La (q + x {1) -k) X s a }j a 
= Jd'qeM-Dq'-F^k^yKd + z^xs^ 

[-2q-ia 2 xk}[y La (q-x^k) X S a }j a (22) 

Altough it is possible to pursue the calculations using (^), experimentally all 
the known transitions exhibit either L a =0 or Lf,=0. The resulting formulae look 



11 



much more sympathetic in those cases and we report here only these special cases 
since only them will be applied in the next chapter. One important difference, 
as compared to the LWLA, is that both the electric and magnetic terms of the 
operator contributes to a given transition. Let us present the result for L b = 0. 

(M A ^ BLb=0 )^ = S Sb ,j b [£^(A; x) + M^(A; x)} (23) 
The term E^A^x) comes from the electric part of the operator and writes: 

S^A- x) = 6 Sa , Sb ^|^4r(5/2) < L a , M a - M b] J b , M b \ J a , M a > (24) 

< 1, -/i; L a - 1, M a - M b + fi | L a , M a - M b > {-lYy La -i,M a -M b +^) 

' < ex > (x( 1 )) L °- 1 exp(-f( 1 )fc 2 ) La < e 2 > exp(-f( 2 )fc 2 ) ' 

mi D 5 / 2 +( ' m 2 D 5 / 2 



Note that this term vanishes if L a = 0. 

The term Ai^{A; x) comes from the magnetic part of the operator and writes: 
M,(A;x) = (-l) 1+s "2^r(3/2)5 a I * ^ X { 2 \ (25) 



1/2 S a J ; 



b 



< ei >xW L «eM-F^k 2 ) ]La +s a +j b <^> x (2) L ° exp(- F^k 2 

mi D 3 /2 + ( ) " " m2 D3/2 

< LaVaSaVa \ J a M a >< S a (J a lv \ J b M b X \v\o | 1/i > y la (k)^ LaMa (k) 



The case L a = looks very similar, but one has to be very careful with the 
phases. In this case the electric part is given by 

E^B- z) = {-I^-^EIJA = B-x = z) (26) 



meaning that in expression (|24|), one has to change all quantities relative to A 
by the corresponding ones relative to B, change the sign of /i, change x by z, 
take the complex conjugate and multiply by a given phase. The expression for 



Aifj,(B;z) is obtained with the same prescription as (26) 



If one admits more than one gaussian function in the expansion of the wave 



function the quantities defined in (21) depends on which terms are retained and 
must be written more explicitly i.e Dij = (Ai + Bj)/2. To obtain the complete 
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expression corresponding to the exact wave function, one must take care of this; 
for example one must make in ( p4| - p5| ) the following replacement: 

( X (i))^-i exp (_^(i) fc2) " AB ( x ^)^ eX p(-F^) 

£)5/2 C * C J" 5/2 ^'^ 

and similar replacements everywhere. 



4 The phase space 
4.1 density of states 

The density of states is obtained with periodic conditions in the quantification 
box. The treatment can be found in any textbook. The density of states by 
energy unit and by solid angle unit is given by 

VE 2 

^ = wr (28) 

One has to calculate the matrix element < Bj\Hi\A > 2 from (|Io|). The best 
way is to use spherical components for the vectors and Racah algebra to deal 
with the corresponding expressions. One has then to sum over the polarisations 
of the photon and of the final meson and to average over the polarisations of 
the decaying meson. In order to simplify the notations, let us introduce the 
quantity X(E) by: 

X ( E )h b ,-k = ^Y, E \<Bi\H t \A>\ 2 (29) 

J A \=±XM a ,M b 

The decay width is given by the golden rule: 

r = J dE 6(E f - Ei) J dn 2ir p{E,tt)X{E) (30) 

In the litterature, one finds different formulae for the width depending upon 
how is treated the energy Dirac term in ( |30| ) . This is known as the phase space 
factor $. Explicitly, one writes : 

T = ^(E )VX(E ) (31) 



As this must be, the quantification volume V, appearing in (31), cancels with 
the one present in X(E) as seen in (|Io|). Eq denotes the energy value fulfilling 
the equation Ef = Ei. 
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4.2 relativistic phase space 

In this case the energies Et and E^ are given by their relativistic expressions. 
Taking into account the fact that the momentum of the final meson is opposite 
to the photon momentum, the Dirac factor is S(Eb + E — m a ) = 5{\Jml + E 2 + 
E—m a ). The integral in (|3^) is performed with the usual rules on delta functions 
to give : 

, s E%E B {E ) m 2 - ml 

$(Eo) = -g- BK 0) ; E = -i b - (32) 

7r m a 2m a 

4.3 Non relativistic phase space 

In a non relativistic treatment, the energies are related to their momenta by the 
classical expressions. Alternatively, one can make the approximation Eq << 
m a , mh in the relativistic phase (p2|). 



$(E ) = ^ ; E =m a -m b (33) 

TT 

Thus, as compared to the relativistic expression, the non relativistic phase space 
differs by two effects. The energy of the photon is equal to the energy difference 
between the resonances (as in nuclear physics) and the term Eb J 1 E "* > is equal to 
unity. 

4.4 mixed phase space 



The mixed phase space is sometimes used in litterature [13 . It consists in 
retaining the relativistic value (|32] ) for the energy Eq - because the effect should 
be important in the meson sector - but the non relativistic value ( p3| ) for Q{Eq). 
Thus the mixed phase space is based on : 

HEo) = ^ ; E a = V Ar^ (34) 
tt 2m a 

5 Total widths 

5.1 long wave length approximation 

With expression (|l5|) for the amplitude and ([n]) and |3^) for the phase space, 
the total width for electric transition is obtained after some calculations: 



f A—'B'y = d S a ,S b 4-Eo ( J 

3 m a mi m2 
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L a J b < L a 010 | L b > 



dpR niLi (p)R naLc 




(35) 



In the same way the magnetic transition width results from the expression 
for the amplitude. It looks like : 



1 A-> 



= SL a ,L h 2aE; 



-2^2-2 
S a Sb Jb 



E B {E ) < ei > 



/ dpR nb L b (p)R 

Jo 




(36) 



These formulae hold for a relativistic phase space. The modification due to 
the use of a non relativistic or a mixed phase space results obviously from the 
discussion of the previous chapter. Simplified expressions applied to specific 
transitions seen experimentally are relegated to the appendix. 
Very often, the electric transition is calculated using the dipole approximation. 
The advantage is that it is not necessary to have the meson wave function 
in momentum representation but in coordinate representation which is a more 
natural scheme, while avoiding a complicated nabla operator. However one must 
be very careful in handling it because it relies on some approximations which are 
not always justified. Indeed the dynamical factor appearing in the amplitude 
is proportionnal to : (B \ p \ A). The widely used trick is to remark that p is 
itself proportionnal to [ H m t, r ] so that (B | p | A) = i/Lt(m a — rrib)(B | r | A), 
where \x is the reduced mass of the meson system (this is sometimes known as 
the Siegert theorem). 

Thus the only change to the previous derivation of the decay width for the 
eletric term in equation ( |35| ) is the replacement of J °° dp Rn b L b (p)Rn a L a {p)p 3 
by /i(m a — nib) f dr Rn b L b (r)Rn a L a (r)r 3 . But, in so doing, one must be sure 
of two things : i) the quark-antiquark potential does not depend on velocity 
(this is the case in our calculations) ii) the meson wave functions are the true 
eigenstates of Hi nt . Even in this case (m a — rrib) is the calculated value with a 
given potential and not the experimental value. 

Moreover, we find also a very popular additional approximation which con- 
sists in the replacement of Eo(m a — mb) 2 appearing in the width by Eq. This is 
acceptable in a non relativistic situation but not in a relativistic one. All these 
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approximations are of course the remnants of the theory applied in atomic or 
nuclear physics. We also discuss this point in the next section. 

In the rest of the text this approximation will be called the dipole approxi- 
mation (DA). 

5.2 general case 

We now come to the expression of the width in the general case, but with the 
wave functions expanded on gaussian terms, as discussed in details previously. 
We present the results only in the case Lb = 0. The results for L a = are easily 
obtained from these ones with the correct replacement (|2^) and the modification 
due to spin average. We do not want to enter into too much details because the 
calculations are essentially a tricky application of Racah algebra. We report the 
result below : 

T A ^ Bl = 5 Sh , Jb 2a EB ^ Eo) [££ + £M + MM] (37) 
m a 

The terms ££,£M, MM come from the electric-electric, electric- magnetic, 
magnetic-magnetic part in the square of the amplitude. One sees that a transi- 
tion is no longer of purely electric or magnetic type, but a mixing of both with 
interference effects. They are given by : 

£M = ^, gb (-i)^+ 5 ° +j °^ +i r(3/ ^ (5/2 V ^(^ + i)g;i:aJ b 

J 1 1/2 1/2 { J L a S a J a \ fc 2L a + l^ ( g g) 



1/2 S a J J I Jb L a 1 

mm = Inm's'L.'j' | J 2 ^ ^ \ ki^v 



2 



(39) 



where the dynamical factors A and B result from the meson wave functions : 

A = < 61 > y c ^(4') £j ' lex p(-4 )fc o) 

mi ^ n 5 / 2 

hJ ^ij 

1,3 ij 



16 



B = 



+ (-1) 



< ei > 



L a +S a + J b 




(41) 



Let us remark that if L a = (transition from S state to S state), the terms 
££ and £M vanishes, and the transition is purely magnetic. Now, we have in 
hand all the tools to perform exact calculations. The formulae look complicated 
but they simplify a lot for transitions of experimental interest. Such simplified 
forms are presented in the appendix. 



We want to stress at the very beginning that, due to the fact that the wave 
functions have been determined with given potentials and that the transition 
operator is perfectly well defined, all the results presented in this section (except 
the very last subsection concerning mixing angles) are free from any adjustable 
parameter. They are thus a very good tool for exploring in detail the drawbacks 
of the formalism. 

6.1 Some comments on the mesons 

In the various potentials under consideration, we do not consider instanton con- 
tribution which can mix strange and ordinary sectors. So, in order to describe 
correctly the scalar 77 mesons, we must introduce by hand some flavor mixing 
angle. Here we consider the idealized case of equal mixing between n flavor 
(isospin doublet u or d) and strange flavor s, which gives: 



Some other comments are in order. The QED conserves the flavor of the 
particles at the vertex (the radiative transitions with flavor change b — ► sj are 
intensively studied recently see [wj but they need penguin diagrams that we 
do not consider here); this implies that the quark content of meson B is the same 
as the one in the initial meson A (in our elementary decay process exhibited in 



6 The results 



v = 




' (nn)i =0 + ss 
With those presciptions the r] and rj mesons are orthogonal. 
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figure (En). Nevertheless the experimental data show a non zero decay width for 
the following transitions concerning vector mesons: cf> — + lo^i and — > pj. For a 
meson taken to be a pure ss, as usually prescribed, this seems to indicate that 
flavor conservation is violated. Instanton effects cannot be advocated because 
they do not play any role for a spin triplet. Those reactions result from more 
complicated processes. One can imagine for instance an annihilation of the qq 
pair in meson A into one or several virtual gluons and creation of a new qq pair 
in meson B (of possibly different flavor) that interacts and emits a real photon. 
This is possible only for neutral flavor mesons. Indeed such transitions (with 
change of flavor) does not occur in kaons for example. 

A possible way to take into account phenomenologically those kinds of tran- 
sitons, in our elementary process, is to include in the wave functions of the to 
a strange flavor part, or/and in the </> wave function a n flavor part. A diffi- 
culty immediatly appears in the case of p resonance since it is isovector while a 
strange flavor can only create an isoscalar. For the 10 and cj) one could introduce 
a mixing flavor angle 9f. The physical mesons are now a combination of the 
ideal (f>o = ss and lo = (nh)j =0 . 



Now with this modification, the transition <fi — > u>^f can be understood as 
the result of two contributions : 4>{ss) — > uj(ss)"f and cj)(nn) — » Lo(nn)j. For the 
4> — > p"f only the n flavor part of the <j> contributes. 

From the electromagnetic point of view, only the charge is conserved, that 
is the projection of the isospin but not the isospin value. We can imagine an 
isospin mixing angle 9i between neutral p and lo. As a consequence p couples 
to (j> at the second order. Another evidence for this possible mixing angle is 
the discrepancy between the neutral and charged channel of the p — » n^f. This 
isospin mixing could be understood by the near degeneracy of p and lo masses 



Our mixing angles are chosen in order to recover the usual prescription for 
9 = 0; this does not correspond always to what can be found elsewhere. For a 
great part of our study, we work in the usual scheme and do not consider mixing 
angles. We will drop those presciptions at the very end. 
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6.2 Influence of the quality for the wave function 

We investigate the quality of our results as the number of gaussians terms N in 
the expansion (^) of the wave function varies from N = 1 (approximation) to 
N = 5 (exact wave function). For this study, we use a relativistic phase space, 
the AL1 potential and focus on the transition p + — ► 7r + 7. The convergence 
properties on the energies have already been discussed in Table ^ . The two 
formalisms - LWLA and the general case - are employed. The results are shown 
m the table (§). 

Table 3: Variation of the decay width (in keV) as function of the number of 
gaussian terms used in the description of the meson wave functions. Long Wave 
Length Approximation (LWLA) and general case (Beyond LWLA), relativistic 
phase space and AL1 potential are reported. 



p+ —> 7T + 7 


1G 


2G 


3G 


4G 


5G 


exp. 


LWLA 


50.98 


56.57 


56.98 


56.98 


57.01 


67.82 ± 7.55 


Beyond LWLA 


44.54 


48.25 


48.50 


48.47 


48.51 


67.82 ± 7.55 



First of all, one can see a sizeable improvement going from N = 1 to N = 2 
for both studies. Limiting oneself to N = 1, as is very often done in literature 
can lead to 10 % error. For N = 2 to N = 3 there is still a further small 
variation. For N = 3 the convergence is achieved and a stable result is reached; 
no further improvement is got passing from N = 3 to N — 5. This is consistent 
with figure ([!]). This means that 3 gaussian terms are enough to describe the 
wave function in a correct way. Nevertheless, for the rest of our study, every 
calculation is performed with N = 5 and the result is considered as the exact one 
concerning the wave function. This is particularly important in the case of radial 
excitations (with one or several nodes in the radial part), since obviously N = 1 
cannot explain such a state and even N = 2 could be a crude approximation. 
Moreover, since our results are analytical, computation with N = 5 is practically 
as fast as the N — 3 case. 

We would like to point out that even if we used only one transition ( 3 S' 1 — > 
^^Sq) to test the quality of the wave function, the same conclusions are valid 
whatever the transition under consideration. We know by experience that N = 
5 is always a very convenient tool to describe the exact wave function. The 
fact that the above transition is better reproduced in LWLA than in the exact 
formalism is not general and is commented later on. 
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6.3 The phase space 

To study which phase space, among the three ones proposed in the last section, 
is the most appropriate for those calculations, we used the AL1 potential, exact 
wave functions and LWLA. In any case, we showed in other works that all 
phase space considerations must be done with the experimental values of the 
resonance energies. The transitions 3 Pj — > 1 5'o and 1 P\ — > 3 Si that are seen 
experimentally cannot been explained in the LWLA and are not reported in this 
part. 

The results for the relativistic phase space (RPS), the mixed one (MPS) 
and the non relativistic one (NRPS) are compiled in the Tables A global 

glance at the quality of the results vs the phase space immcdiatly reveals that 
(NRPS) is not a correct prescription, the theoretical results being very often too 
large, except may be for heavy mesons for which a non relativistic treatment is 
acceptable. The problem comes essentially from a wrong determination of the 
momentum fco = Eo, that enters both in the J^- factor (here taken as unity) 
and in the dynamical amplitudes. There exist some important discrepancies 
between the relativistic and the non relativistic momenta (this could be seen 
in the second and third columns of the different tables). We remind you that 
the experimental values of the meson masses are used. The MPS supports this 
conclusion; it uses the relativistic impulsion but a ^j- factor equal to one (as 
in the non relativistic phase space) and the results are much better than NRPS 
ones, and little poorer than RPS ones. So the predicted decays widths in NRPS 
are essentially spoiled by wrong values of momentum. The change in the results 
going from MPS to RPS is proportional to the factor So looking at the 

fourth column of the tables, one sees that the biggest deviation is around a 
factor 2 and concerns the transitions : p + — > 7r + 7 or ui — > 7r°7. The use of 
RPS is quite satisfactory for the whole bulk of data, indicating that the wave 
functions are not completely crazy. 

In the fourth column we also see that -§2. = 1 is a justified approximation 
for the heavy mesons, this is due to the fact that the mass spectrum is here 
more dense than in the light sector. 

With those values of momenta and considering their important values for 
some transitions, it is also interesting to be aware that the LWLA could not be 
always a justified formalism. In the case of important mass differences between 
initial and final state, the recoil term could not be omitted. 

In view of this discussion, the rest of our study will be done using a relativistic 
phase space. 
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6.4 Dipole approximation 

In this part we want to discuss the dipole approximation (DA) for electric tran- 
sition. A relativistic phase space is used, but we study two different prescrip- 
tions suggested previously. DAI is the expression resulting from the formalism 
: it appears a term Eo(m a — nit,) 2 in phase space; DA2 is obtained with this 
term replaced by Eq. The results are presented in the table (|^). It is diffi- 
cult to draw some conclusions , mainly because of lack of precise experimental 
data. We first remark a small difference (less than 15%) between the DAI and 
DA2. This is due to the fact that almost all exploitable data concern heavy 
mesons for which the two approximations tend to be the same. For example, 
the transition Xb2(2P) —> T(lS)j gives a momentum value of E — 776 MeV 
for m a — mb = 808 MeV, that is a 8% variation. Indeed only two transitions 
deal with ordinary quarks : /i(1285) — * ( o°7 and &i(1235) + — > 7T+7 and in this 
case the difference is appreciable. 

Generally speaking the DA gives better agreement than LWLA on weak 
transitions (around 1 keV) but this is not significant. On the contrary the 
results are worse in the case of strong transitions. Although the DA is very 
popular in atomic and nuclear physics - a domain where it can be applied safely 
- it should be avoided in the meson sector. 



6.5 General study for three different quark-antiquark po- 
tentials 

In this part we present the most sophisticated calculations in this framework. 
We go beyond LWLA, use relativistic phase space and exact wave functions. The 
objective is twofold. First to test under which conditions a general treatment 
is necessary as compared to the often used LWLA, and how big could be the 
difference. Second, by testing the three quark-antiquark potentials proposed 
in section 2A, to see whether the results are very sensitive to the dynamics of 
quarks inside a meson. The results are presented in tables (0§]9|)- 

Let us compare this general case with the LWLA. Several comments are in 
order. In a number of cases there is no electric-magnetic mixing. For L a = 
Lb = only the magnetic term remains while a pure electric term remains for 
1 Pi — > 1 So- In this case the LWLA gives always a larger value; this can be shown 
theoretically if the wave function has only one gaussian component. Since in 
general one component is dominant it is not surprising that this property persists 
even in more realistic situations. Although this is not always the case, LWLA 
often gives better (as compared to experiment) results; this means that either 
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the wave functions are not so good or that something is still missing in the 
theory. 

The transitions 1 P i — > 3 Si and 3 Pj — > 1 S'o are also purely magnetic but 
are completely forbidden in LWLA. Our complete treatment predicts them with 
the right order of magnitude. 

Lastly 3 Pj — > 3 Si and 3 Si —> 3 Pj transitions show interference between 
electric and magnetic contributions while they are purely electric in LWLA. 
Since the electric part remains largely dominant there is no big difference be- 
tween both formalisms except in the case of /i(1285) — > pj which is greatly 
enhanced in the right direction. 

An overall look at those tables shows us the similarity of the decay widths 
resulting from the AL1, API and Bhaduri potentials. The results obtained with 
AL1 generally lie between those of Bhaduri and API. The predicted values 
coming from the AL1 potential are smaller than the API ones. This could be 
related to the different asymptotic behavior of the potential at long range. The 
confinement being slower in the API potential, the spread of the wave function 
is more important and contributes more in the spacial integration. Globally no 
potential is really more suited than the other for those calculations, although 
Badhuri's one seems a little poorer. The trends are essentially the same and 
when one potential gives too low (or too high) value, so do the others. The 
agreement with experiment is satisfactory for all types of transitions, giving 
indication that we are in the good track for the description of mesons. The 
discrepancy with experimental situation exceeds 50% very scarcely and this can 
be considered as encouraging but not completely satisfactory. 

The fact that three different wave functions give more or less the same trends 
(although there can exist 20% differences) shows that the quality of the wave 
function is not responsible entirely for the discrepancy with experiment. More- 
over LWLA as an approximation should give poorer agreement than a complete 
exact treatment. This is not the case (except of course when it gives a null 
result). This proves that something is still missing in the formalism although 
the present day calculations provide the dominant contributions. 

Now let us have a closer look to some interesting transitions. A special one 
concerns the neutral and charged decays of the p into the ttj. Experimentally 
the decay width for the charged channel is 68 keV whereas the measure for 
the neutral is 102 keV. In our calculation the small difference between those 
two channels comes only from the tiny difference between the experimental 
masses of the tt + and the 7r°. The decay width has the same expression for 
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those transitions due to the term 



mi ni2 

meson composed of a single flavor. That is [| — i] 2 = i for the charged channel 



in expression (|25| ) appearing for 

1-12 _ 1 
3 3 J 9 



and [i — -jf] 2 — g for the neutral one. So where does this important variation 
between those two channels come from ? First we have to point out that given 
the large uncertainties 67.82 ± 7.55 and 102.48 ± 25.69 the two values are nearly 
compatible with 76 MeV. This means that there may be no problem with those 
channels except an experimental one ! Nevertheless if we rely more deeply on the 
experimental values a possible insight could come from the u> — > ir°j transition, 
which is identical to p° — > 7r°7 but with an enhancement by a factor 9 due 
to isospin. The experimental value 715 keV is roughly in agreement with this 
point. So even a small isospin mixing between the u and the p° could increase 
sufficiently the decay width to explain the data. This hypothesis will be tested 
in the next part. 

We remark the same variation for the K* decaying into the K meson and 
for the B* into the B. Yet this time the isospin factor explains this variation. It 

appears a factor (3^— + j^r) 2 f° r the neutral channel and a factor 3^r) 2 

for the charged one in the LWLA (it is more complicated to estimate the ratio 
of the strange with the isospin doublet masses in the general formalism ) . Using 
the experimental data and making the approximation that the matrix elements 
in those two channels are identical (except the isospin dependence) , we find the 
relation : m s = 1.24to„. In our potentials this ratio is 1.83, 1.78, 2.00 for 
the AL1, API and Bhaduri respectively. 

Concerning the transition a? — ► 7r°7, the decay width is zero; this is due 
to the fact that for the 3 Pj — > 1 So composed of a single flavor the width is 
propotional to + = 0. 

In the potentials used, there is no isospin dependence so the p, uj have the same 
radial part of the wave function, and the same remark is true for the tx and rj n 
(no instanton effect). This means that for the transitions number one to five of 
the 3 Si — » 1 5o tables, the integral part is identical. Now since the experimental 
hierarchy of the decay widths is reproduced in our calculation this could indi- 
cate that we used a correct prescription; moreover since the numerical value are 
also reproduced this ensures the quality of our wave functions. 
It is not sure that the decay of the D s i(2536)* + into D**y has been observed 
experimentaly but our result for this width tends to prove that it should be seen 
experimentally. 
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6.6 Mixing angles 

If the wave function is composed of two parts as in the 77 mesons (flavor mixing) , 
or in the p (isospin mixing with the to), some precisions are needed in the formal- 
ism. In the case of 77, the wave function can be written: |^ >= \^ Vri (nn; I = 
0) > — l^r/s >• In our study there is no instanton effect so we put by hand 
a mixing of 50% between the two flavors. That is we calculate separately the 
states fj n and rj s both normalised to one and we have |\t^(nfi; I — 0) >= -2)=, 
\^Vs >= A possible difficulty is that those states do not have the same mass 
in order to calculate the phase space; nevertheless as we take the experimental 
value there is no difficulty 

From a general point of view we write >= j^i > ±| v I / 2 > where 1 and 2 
denote the two flavor (or isospin) components of the wave function. We have to 
calculate: 

M(A -> Bi) =< A\M\B >= (< Ai\± < A 2 |)M(|Bi > ±\B 2 >) 

so there could exist 4 components, and therefore some interferences. In the case 
of flavor mixing not all the terms will contribute due to flavor conservation. 



As we have said in section 6.1, the </> could decay into the u>, the p and the w. 
This could be incorporated to our decay process by two ways: an isospin mixing 
(u> and p) or a flavor mixing (</> and lj). This study is done with the general 
formalism and a relativistic phase space. Because of the similarity of the results 
obtained via the three potentials, it is sufficient to perform the calculation with 
only one potential, here the ALL 

6.6.1 Flavor mixing 

In this part, we investigate the mixing between the <fi an d the u>. Technically 



we use the prescription of section 3.L Now we have to find an appropriate 
value of the mixing angle 9 f and for that we rely on experimental data. A good 
candidate is the transition <j> — > ir°j which is possible only through a flavor 
mixing. Only the nn flavor part of the <f> contributes to the decay. We find a 
small mixing of Of = 4.5 degrees. We could have choosen oj — > 777 to determine 
the angle value but the transition including a 77 meson are not very appropriate 
because of its flavor mixing which could generate some interference terms. The 
u — » 7r°7 transition is no more suited for this, even if only one term (u> n n) 
contributes because the value of a pure nn meson is smaller (459.30 keV) than 
the experimental value (714.85 keV) so including a ss part to the u which will 
not contribute to the decay could only decrease the width. 
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The transitions modified in consequence are presented in table ( [10| ) . First of 
all, we see that for the transitions already allowed, there is no superority taking 
into account mixing angles. Some transitions like the lo — > 777, 0(1020) — > 
777 are deteriorated whereas the 77' (958) — > ^7 is improved. In fact the only 
improvement is the possibility to calculate transitions like the /i(1285) — > 07, 
<f> — > fo (980)7, 4 1 ~~ > a Q (980)7 which would be forbidden otherwise. Nevertheless 
the mixing is here to small to reproduce the experimental data. 

6.6.2 Isospin mixing 

Here we consider an isospin mixing between the lo and the p°. This mixing 
do not appear for the charged p because it is an Mj = ±1 meson. The angle 
07=8.9 degree is taken to fit the transition p° — > 7r°7, and discriminates between 
charged and neutral transitons. The results are presented in table ([Tl|). 

Considering the poor quality of the results (4 transitions arc deteriorated 
and 2 improved but not reproducing the experimental values), it is clear that 
we are missing something. This could be very well a false angle value, may be 
the p° — > 7r°7 transition results of another process and should not be use to 
determine 9j. 



7 Summary 

This work is a first review of the decay of a meson into another one plus a 
real photon. We have analyzed carefully the different parts of this elementary 
process. First of all we have presented different formalisms; the Long Wave 
Length Approximation, also called the static approximation, where the recoil 
term is neglected. We also investigated the Dipole Approximation, in which the 
spacial integrals are performed in the representation space with the presence 
of the factor (tua — ™_b); another step in this non relativistic approximation 
is the replacement of (rriA -wb) by the photon energy Eg. And finally we go 
beyong the Long Wave Length Approximation. The different formulae for the 
decay widths are presented in the appendix. In view of the results no formal- 
ism seems more suited than the others to describe the radiative decay, all of 
them giving the right order of magnitude, nevertheless the superiority of the 
general formalism comes from the fact that it allows the calculation of electric- 
magnetic interference terms and forbidden transitions in the LWLA and DA 
such as 1 Sq — ► 3 Pj. 

Secondly, we study three phase spaces: relativistic, non relativistic and mixed. 
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This last one uses a relativistic momentum but a non relativistic value of the 
ratio which allows us to show that the important point is a good description 
of the momentum. Anyhow we showed that a relativistic phase space is always 
better and should be used in any circumstance. 

Then we checked the importance of the wave function through the use of three 
potentials: AL1, API and BD. Those calculations do not allow to show the su- 
periority of one of them, the predicted values being of the same quality. For our 
results to be analytical we expand our wave functions as a sum of N gaussian 
terms. We showed that N=3 is sufficient to obtain a convergence of the results 
but in order to be sure to treat the exact wave function, we used everywhere in 
our calculation N=5. We also showed that using exact wave functions is always 
preferable as the one gaussian approximation that is sometimes of common use. 
The difference could be as large as a 20% effect. 

Finally we incorporated some mixing angles in order to calculate some otherwise 
badly reproduced or even forbidden transitions such as <f) ~ > ?r 7. Those angles 
are of two kinds: isospin mixing between the p and u mesons and flavor mixing 
between the (f> and u) mesons. The results are not satisfactory. Except for a 
possible explanation of the important difference between neutral and charged 
decay of the p into pion, the isospin mixing deteriorate the quality of the pre- 
dicted value whereas for the flavor mixing there is no important improvements. 

Although a completely rigourous formalism gives an overall satisfactory 
agreement with experimental data, we gave arguments that, in the framework 
that we considered (NRQM and NR expression for the transition operator) some 
physics is still absent. We think that one can explore two different directions : 
inclusion of relativistic effects both in the wave functions and in the transition 
operator, and introduction of form factors at the quark-photon vertex. Such 
works are under study. 

Appendix 

Long Wave Length Approximation 

In this part, we give the simplified expressions of ( j35| , j36| ) for the transitions of 
experimental interest, using wave functions expanded on gausssian terms as in 
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r( 3 5i - 3 Pj) 
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We see that for the last formula, this transition is null in the case of a meson 
composed of only one flavor due to isospin term. 

Forbidden transitions 

The following transitions are not allowed in the LWLA. 

• tto - 3 P.J 

• 3 Pj - % 
. -> 3 ^i 



General Case 

Here we give the simplihed expression of the total exact width 



27 



L b = 

We adopt the following notation for the decay width formulae: 
T A ^ Bl = S Sb ,j b 2a EB{Eo) [££ + SM+ MM] 

^A-^B-y = 

££ = 0; £M = 

32 11 mj ^ D l m 2 ^ 



3 Si - % 

= 0; £7W = 

(A B -_F (1) B 2 
mi ^ n5 m 2 4V 

3 5i - 3 Si 

££ = 0; £7W = 



(1) 2 

MM = ^ 3 [<£i>y ^g^ + <^>y;«. 

48 ' 1 mi ^ r>§ m 2 4V 

«J u i,j hJ 



• ^1 - % 

£M = 0; MM = 
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cc = I_ E / < ex > v cfcfe-^^ <e2> cfcfe-^?^ * 



f £ = 0; £M = 



= 0; £M = 

- J = MM = 



- J = l MM = ££ 5 ^Eij™4 + ^£^ 



- j = 2 mm ^ (^E i; ,™4 + ^£ 4J ™- F 



• 3 Pj - 3 5i 

- J=0 



= — £ -4 2 (£o); £M = £ 3 *4(£o)£(£o); MM - ^ i^ 2 (£o) 

- J=l 

55 = ^ £; ^ 2 (^o); £M = ^ ^(i^W); MM - ^ E 5 Q B 2 (E Q ) 



J=2 



f 5 = ^ £ -4 2 (£o); 5M = E 3 Q A(E )B(E ); MM = ^ EqB 2 (E ) 



with: 
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m i n= m 2 ^ n = 

A B (1) -F (1) E 2 A B (2) -F (2) E 2 

atjP^ <ex> ^cfcfxl/e » ° <e2> cfcfxl/e » 

b(h ) = > 3 > 3 

L a = 

We adopt the following notation for the decay width formulae: 
r A ^B 7 - <5 Sa ,./ a 2a^^- [££ + £M + MM] 

• % - ^1 

£M = 0; MM = 

(1) 2 (2) 2 ' 

rr 3^ _ / < ei > ^ cf e" F « £ » < e2 > _ cfcfe^ B « 
32 \ m i ft 7712 Tl D 2 



5 £ = 0; £M = 

J = jMM = 

J = l MA4 = ff^ f-^-V / ■ \\ " " +^E^ 



_4 B (1) -F^eZ 

cycfzY, ,'e »J 



D 2 - 

t-,3 



- J = 2 MM = |f^ -^-V + ^E,/ 1 

\ > j 



= 0; £M = 

A B (1) - 



mm = — ei I ^ C1 ^ y ? j Zii ! + <62> y CiC ^ e 



32 \ mi ^ n 2 m 2 z - / n 
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3 S1 - 3 PJ 

- 3=0 

- 3=1 



EE = — E A 2 (E ); £M = — EqA(E )B(Eq); MM = — EqB (E a ) 



3=2 



£ £ = — E A 2 (E ); SM = — EqA(E )B(Eq); MM = — E 5 B 2 (E Q ) 



tlT iTT i IT 

££ = —E Q A 2 (E Q ); £M = -—E$A(E )B(E ); MM = —EIB 2 (E ) 



with: 



-£c?e~ F *P E ° <r Pn -> c A c¥e~ F ^ E ° 



A(E ) = ^^Y CiC i e „" "_<*>Y**±_ 

m-t t^o rrio ^ 



b{e ) = 1 1 3 lJ - <e2> 1 a zii 

1,3 1,3 
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Figure 1: tt and p radial wave functions ^|=r f° r N=l to 5 gaussian terms 
compared to the exact ones. 
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Pa - k 
~Pa = 



mi+m-2 
mi+m-2 



k + Pb 
- Pb 



Pa 
~Pa 



Pa 
~Pa ■ 



mi +m 2 

k ^ 



k + Pb 
k - Pb 



mi +m 2 



Figure 2: Schematic representation of the elementary decay process with mo- 
mentum values from the rest frame of the inital meson. (1) denotes the emission 
of the real photon by the quark and gives rise to AfW and (2) by the antiquark 
and gives rise to M^> . 
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Table 4: Decay widths obtained through the Long Wave Length Approxima- 
tion, relativistic (rel), mixed and non relativistic (nrel) phase spaces, AL1 po- 
tential. The momentum fco = Eq is expressed in MeV whereas the decay widths 
are in keV. The mixed phase space uses a relativistic momentum contrary to 
the non relativistic one which also uses a ratio ^f- equal to the unity. 



3 Si - % 


transition 


jLreZ 


unrel 
K 


Eb_ 

m A 


rel. 


mixed 


nrel. 


exp. 


p + — > 7T + 7 


372 


630 


0.52 


56.98 


110.32 


535.23 


67.82 ± 7.55 


p° -» 


373 


634 


0.52 


57.24 


111.05 


547.02 


102.48± 25.69 


p° -> 777 


190 


222 


0.75 


49.79 


66.12 


105.52 


36.18 ± 13.57 


W — ► 7T°7 


380 


648 


0.51 


543.44 


1055.48 


5238.71 


714.85 ± 42.74 


w — > 777 


200 


235 


0.74 


6.37 


8.56 


13.96 


5.47 ± 0.84 


</>(1020) -> 777 


363 


472 


0.64 


45.90 


71.26 


157.04 


55.82 ± 2.73 


^(1020) -» 7/(958)7 


60 


62 


0.94 


0.30 


0.32 


0.35 


0.53 ± 0.31 


if* (892)° -> K°7 


310 


398 


0.65 


111.19 


169.95 


361.34 


116.15 ± 10.19 


A"* (892)+ -» if+7 


309 


398 


0.65 


97.47 


149.21 


318.28 


50.29 ± 4.66 


L>*(2007)° -> L»°7 


137 


142 


0.93 


35.08 


37.65 


41.95 


< 800.10 ± 60.90 


L>*(2010)+ -» L»+7 


136 


141 


0.93 


2.68 


2.87 


3.19 


^ 1 AA + 2 - 7 5 




139 


144 


0.93 


0.28 


0.29 


0.33 


< 1789.80 ± 47.50 


B*+ -» B+7 


46 


46 


0.99 


0.97 


0.98 


0.99 


seen 


£T° -> B°7 


45 


46 


0.99 


0.28 


0.29 


0.29 


seen 


J/* -» 77 c7 


115 


117 


0.96 


1.85 


1.93 


2.04 


1.13 ± 0.35 


^(25) - r? c (15)7 


639 


706 


0.83 


2.15 


2.60 


3.52 


0.78 ± 0.19 




transition 


urel 


unrel 


E B 
m A 


rel. 


mixed 


nrel. 


exp. 


r/(958) -» /9°7 


170 


188 


0.82 


116.80 


141.99 


193.73 


61.31 ± 5.51 


7/(958) -> W7 


159 


175 


0.83 


10.81 


12.97 


17.29 


6.11 ± 0.78 



35 



Table 5: Same as Table (@). 



3 Si -> 3 Pj 


transition 


h,rel 
K Q 


h,nrel 
K Q 


E B 
m A 


rel. 


mixed 


nrel. 


exp. 




261 


271 


0.93 


17.70 


19.04 


19.77 


25.76± 3.81 


1>(2S) Xci(lP)7 


171 


175 


0.95 


35.75 


37.49 


38.40 


24.10 ± 3.49 


?P(2S) -> Xc2 (lP)7 


127 


130 


0.97 


44.90 


46.51 


47.35 


21.61 ± 3.28 


T(25)^ X bo(lP)7 


162 


163 


0.98 


0.43 


0.44 


0.44 


1.89 ± 0.53 


T(2S)- XM (1P) 7 


130 


131 


0.99 


1.04 


1.05 


1.06 


2.95 ± 0.61 


T(2S)^ Xb2 (lP) 7 


110 


111 


0.99 


1.47 


1.48 


1.49 


2.90 ± 0.61 


T(35)^ X bo(2P)7 


122 


123 


0.99 


0.68 


0.69 


0.70 


1.42 ± 0.25 


T(35)^ X61 (2P) 7 


100 


100 


0.99 


1.67 


1.69 


1.69 


2.97 ± 0.43 


T(3S)^ Xb2 (2P) 7 


86 


87 


0.99 


2.42 


2.44 


2.45 


3.00 ± 0.45 


3 Pj -> 3 Si 


transition 


h.rel 
K 


h,nrel 
K Q 


E B 
m A 


rel. 


mixed 


nrel. 


exp. 


/i(1285) -» A 


410 


513 


0.68 


939.73 


1381.81 


1727.12 


1296.00 ± 295.20 


Xco(lP) J/V(15)7 


303 


318 


0.91 


233.40 


256.15 


268.67 


92.40 ± 41.52 


Xd(lP) -► J/V(15)7 


389 


414 


0.89 


292.29 


328.75 


349.33 


240.24 ± 40.73 


x c2 (ip) - JMish 


430 


459 


0.88 


319.01 


362.85 


387.90 


270.00±32.78 


X60 (1P) - T(15) 7 


392 


400 


0.96 


28.82 


30.01 


30.63 


seen 


X6i(lP) - T(15) 7 


423 


432 


0.96 


31.03 


32.42 


33.14 


seen 


X62(1P) - T(15) 7 


442 


452 


0.96 


32.37 


33.88 


34.67 


seen 


Xb0 (2P) - T(1S) 7 


743 


772 


0.93 


12.12 


13.07 


13.58 


seen 


X60 (2P) - T(2S) 7 


207 


209 


0.98 


13.07 


13.34 


13.47 


seen 


X6i(2P) -> T(15) 7 


764 


795 


0.93 


12.45 


13.45 


13.99 


seen 


X61 (2P) ^ T(2S) 7 


229 


232 


0.98 


14.47 


14.80 


14.97 


seen 


X62 (2P)^T(1S) 7 


776 


808 


0.92 


12.63 


13.66 


14.22 


seen 


Xb2 (2P)^T(2S) 7 


242 


245 


0.98 


15.27 


15.63 


15.82 


seen 


1 Pi — > 1 Sq 


transition 


Urel 


h,nrel 


E B 
m A 


rel. 


mixed 


nrel. 


exp. 


61 (1235)+ -> 7T+7 


607 


1090 


0.51 


209.95 


414.55 


744.58 


227.20 ± 58.60 



36 



Tabic 6: Dipole approximation in the LWLA, AL1 potential and relativistic 
phase space. Two approximations are presented depending upon the replace- 
ment (DA2) or not (DAI) of (ras — tua) by Eq, see text. 



3 Si - s Pj 


transition 


fco[MeV] 


E B 
m A 


DAI 


DA2. 


exp.[keV] 


4>(2S) -> Xco(lP)7 


261 


0.93 


49.92 


53.80 


25.76± 3.81 




171 


0.95 


43.43 


45.57 


24.10± 3.49 


^(2S)^ Xc2 (lP)7 


127 


0.97 


30.23 


31.32 


21.61± 3.28 


T(2S)^ X 6o(lP)7 


162 


0.98 


1.52 


1.54 


1.89 ± 0.53 


T(2S)-> X61 (1P) 7 


130 


0.99 


2.34 


2.37 


2.95 ± 0.61 


T(2S)^ Xm (1P) 7 


110 


0.99 


2.39 


2.41 


2.90 ± 0.61 


T(3S) - X6o(2P) 7 


122 


0.99 


1.65 


1.67 


1.42 ± 0.25 


T(3S) ^X6i(2P)7 


100 


0.99 


2.67 


2.69 


2.97 ± 0.43 


T(3S)^x 6 2(2P) 7 


86 


0.99 


2.91 


2.93 


3.00 ± 0.45 


3 P/ - 3 Si 


transition 


k 


E B 
m A 


DAI 


DA2 


exp. 


/i(1285)^p° 7 


410 


0.68 


824.56 


1288.17 


1296.00± 295.20 


Xco(lP) - J/V-(1S)7 


303 


0.91 


144.59 


159.06 


92.40 ± 41.52 


Xd(lP) - J/V(1S) 7 


389 


0.89 


298.24 


336.75 


240.24 ± 40.73 


Xc 2 (lP) - J/V(15) 7 


430 


0.88 


396.53 


453.17 


270.00 ± 32.78 


X6o(lP) - T(1S)7 


392 


0.96 


20.36 


21.21 


seen 


X6i(lP) - T(1S) 7 


423 


0.96 


25.59 


26.74 


seen 


X6 2 (1P) - T(1S) 7 


442 


0.96 


29.14 


30.52 


seen 


X6o(2P) - T(1S) 7 


743 


0.93 


11.32 


12.22 


seen 


X6o(2P) - T(2S) 7 


207 


0.98 


9.37 


9.56 


seen 


X6i(2P) - T(1S) 7 


764 


0.93 


12.30 


13.31 


seen 


X6i(2P) - T(2S) 7 


229 


0.98 


12.77 


13.06 


seen 


X6 2 (2P) - T(1S)7 


776 


0.92 


12.89 


13.97 


seen 


X6 2 (2P) - T(2S) 7 


242 


0.98 


15.04 


15.41 


seen 


x Pi - ^0 


transition 


k 


E B 

m A 


DAI 


DA2 


exp. 


6i(1235) + -> 7T+7 


607 


0.51 


84.33 


272.06 


227.20 ± 58.60 
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Table 7: General case, relativistic phase space, AL1, API and Bhaduri's po- 
tentials. For the values of kl el and ^f- see Tables (|^-||)- The electric, magnetic 
and interference terms refer to the AL1 potential. 



3 Si -» 'So 


transition 


elec. 


interfer. 


magn. 


tot(ALl) 


tot(APl) 


tot(BD) 


exp. 


p + — > 7T + 7 






48.48 


48.48 


60.41 


44.07 


67.82 ± 7.55 


p° -» 7T 7 






48.66 


48.66 


60.64 


44.24 


102.48± 25.69 


P° -> ??7 






47.73 


47.73 


60.63 


42.78 


36.18 ± 13.57 


— » 7T°7 






459.30 


459.30 


571.79 


417.85 


714.85 ± 42.74 


uj — > 777 






6.08 


6.08 


7.72 


5.43 


5.47 ± 0.84 


0(1020) -» 7/7 






41.27 


41.27 


44.12 


39.74 


55.82 ± 2.73 


0(1020) -> 77'(958)7 






0.30 


0.30 


0.32 


0.29 


0.53 ± 0.31 


K*(892)° -» if 7 






98.28 


98.28 


116.41 


91.70 


116.15 ± 10.19 


if* (892)+ -> K +7 






79.07 


79.07 


104.46 


71.50 


50.29 ± 4.66 


L>*(2007)° -> D°7 






33.60 


33.60 


41.74 


30.24 


< 800.10 ± 60.90 


D*(2010)+ -> £)+7 






2.48 


2.48 


3.58 


2.09 


<. 1.44 _ .92 


- ^7 






0.26 


0.26 


0.31 


0.22 


< 1789.80 ± 47.50 


s »+ _^ B + 7 






0.97 


0.97 


1.26 


0.84 


seen 


B* B° 7 






0.28 


0.28 


0.36 


0.25 


seen 


J/* -» t; c 7 






1.85 


1.85 


1.87 


1.79 


1.13 ± 0.35 


1>(2S) -» 77.(15)7 






4.97 


4.97 


6.34 


2.44 


0.78 ± 0.19 




transition 


elec. 


interfer. 


magn. 


tot(ALl) 


tot (API) 


tot(BD) 


exp. 


77' (958) -> p°7 






112.90 


112.90 


143.62 


101.09 


61.31 ± 5.51 


t/' (958) -> W7 






10.50 


10.50 


13.36 


9.39 


6.11 ± 0.78 
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Table 8: Same as Table (0). 



3 Si -» 3 Pj 


transition 


elec. 


interfer. 


magn. 


tot(ALl) 


tot(APl) 


tot(BD) 


exp. 


<P(2S) - X co(lP)7 


18.00 


-4.12 


0.24 


14.12 


14.06 


14.94 


25.76± 3.81 


i>(28) xailPh 


36.02 


-1.82 


0.05 


34.25 


34.23 


35.70 


24.10 ±3.49 


i>(28) Xc 2 (lP)7 


45.09 


1.27 


0.03 


46.39 


46.43 


48.07 


21.61±3.28 


T(2S)^ X 6o(lP)7 


0.43 


-0.02 


0.00 


0.41 


0.54 


0.45 


1.89 ± 0.53 


T(2S)- XM (LP) 7 


1.04 


-0.02 


0.00 


1.02 


1.35 


1.13 


2.95 ± 0.61 


T(2S)- XM (1P) 7 


1.47 


0.02 


0.00 


1.49 


1.97 


1.64 


2.90 ± 0.61 


T(3S)^X6o(2P) 7 


0.68 


-0.02 


0.00 


0.66 


0.73 


0.72 


1.42 ± 0.25 


T(3S)^X6i(2P)7 


1.67 


-0.02 


0.00 


1.65 


1.84 


1.79 


2.97 ± 0.43 


T(35) ^Xb2(2P) 7 


2.42 


0.02 


0.00 


2.44 


2.71 


2.64 


3.00 ± 0.45 


3 Pj -» 3 Si 


transition 


elec. 


interfer. 


magn. 


tot(ALl) 


tot(APl) 


tot(BD) 


exp. 


/i(1285)^p° 7 


688.97 


417.42 


126.45 


1232.83 


1376.96 


1209.32 


1296.00 ±295.20 


Xco(lP) - J/V(1S')7 


225.25 


29.20 


0.95 


255.40 


260.24 


271.53 


92.40 ±41.52 


Xci(LP) - JM1S) 7 


275.69 


29.38 


1.57 


306.63 


312.43 


327.31 


240.24 ±40.73 


Xca(LP) - J/V(1S) 7 


297.08 


-38.52 


3.50 


262.05 


266.99 


286.79 


270.00 ± 32.78 


Xw(lP) - T(15) 7 


28.29 


1.78 


0.03 


30.10 


30.85 


32.24 


seen 


X u(lP) -» T(1S) 7 


30.37 


1.11 


0.02 


31.51 


32.26 


33.82 


seen 


X62(1P) T(15) 7 


31.62 


-1.26 


0.04 


30.39 


31.03 


32.82 


seen 


Xw ,(2P) -► T(1S) 7 


12.22 


1.73 


0.06 


14.01 


11.80 


14.69 


seen 


XW ,(2P) -» T(2S) 7 


12.88 


0.43 


0.00 


13.31 


13.52 


14.28 


seen 


X6i(2P) -» T(15) 7 


12.55 


0.94 


0.04 


13.53 


11.38 


14.25 


seen 


X6i(2P) - T(2S) 7 


14.21 


0.29 


0.00 


14.51 


14.72 


15.58 


seen 


Xi, 2 (2P) -+ T(1S) 7 


12.74 


-0.99 


0.05 


11.80 


9.90 


12.60 


seen 


X6 2 (2P) -^T(2S) 7 


14.97 


-0.34 


0.01 


14.63 


14.82 


15.77 


seen 
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Table 9: Same as Table (g). 











'Pi- 


> 'So 








transition 


elec. 


interfer. 


magn. 


tot(ALl) 


tot(APl) 


tot(BD) 


exp. 


61 (1235)+ - 


> 7T + 7 


148.68 






148.68 


152.76 


142.14 


227.20 ± 58.60 


3 Pj 'S 


transition 


elec. 


interfer. 


magn. 


tot(ALl) 


tot(APl) 


tot(BD) 


exp. 


ai(1260)+ - 


+ 7T + 7 






179.53 


179.53 


229.90 


151.63 


seen 


01 (1260)° - 


A 














seen 


a 2 (1320) + - 


■> 7T + 7 






142.01 


142.01 


179.27 


120.87 


299.60 ± 65.71 


'Pi -» 3 Si 


transition 


elec. 


interfer. 


magn. 


tot(ALl) 


tot(APl) 


tot(BD) 


exp. 


D s i(2536)* + 








10.97 


10.97 


11.99 


9.15 


probably seen 
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Table 10: Decay widths obtained with a mixing angle of Of = 4.479° between 
the and u> mesons. Beyond the long wave length approximation, AL1 poten- 
tial, relativistic phase space. The mixing angle is obtained from the — > -K°-f. 





3 Sj 


- x 5 









transition 


fc 


m/t 


theo. 


exp. 


UJ — > 7T°7 


380 


0.51 


456.49 


714.85 ± 42.74 


u> — ► ?77 


200 


0.74 


7.22 


5.47 ± 0.84 


0(1020) -> 7r°7 


501 


0.51 


fitted 


5.617 ± 0.447 


0(1020) -> 777 


363 


0.64 


35.88 


55.82±2.73 


0(1020) -» 7/ (958)7 


60 


0.94 


0.34 


0.53 ± 0.31 




% - 3 S 


1 






transition 




E B 


theo. 


exp. 


r/(958) -» ^7 


159 


0.83 


8.58 


6.11 ± 0.78 




transition 




E B 
m A 


theo. 


exp. 


A (1285) - 07 


236 


0.82 


0.48 


17.76±6.30 




3 Si 


- 3 P., 






0^/o (980)7 


39 


0.96 


0.03 


1.52±0.18 


-> a (980)7 


34 


0.97 


0.23 


< 99 99 
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Table 11: Decay widths obtained with a mixing angle of 6 — 8.88° between the 
p and u> mesons. Beyond the long wave length approximation, AL1 potential, 
relativistic phase space, the determination of the mixing angle is based on the 
transition p° — > n a -f. 









3 Si- 






transition 


k 


m 4 


theo. 


exp. 


p° - 7T°7 


373 


0.52 


fitted 


102.48± 25.69 


p° -> 777 


190 


0.75 


51.57 


36.18± 13.57 


W — > 7T°7 


380 


0.51 


402.85 


714.85 ± 42.74 


w — > 777 


200 


0.74 


1.68 


5.47 ± 0.84 










3 s 1 




transition 


k 




theo. 


exp. 


r/(958) -» 


A 


170 


0.82 


121.99 


61.31 ± 5.51 


r/(958) -» 


CJ7 


159 


0.83 


2.89 


6.11 ± 0.78 














transition 


fc 


Eb 
m A 


theo. 


exp. 


/1 (1285)- 


-P°7 


410 


0.68 


1332.09 


1296.00 ± 295.20 
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